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In this paper we present a class of exact inhomogeneous solutions to Einstein's equations for higher 
dimensional Szekeres metric with perfect fluid and a cosmological constant. We also show particular 
solutions depending on the choices of various parameters involved and for dust case. Finally, we 
examine the asymptotic behaviour of some of these solutions. 
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Usually, cosmological solutions to Einstein's field equations are obtained by imposing 
symmetries [1] on the space-time. One of the reasonable assumptions (in an average sense) 
is the spatial homogeneity. But when we consider cosmological phenomena over galactic 
scale or in smaller scale (detailed structure of the black body radiation) then we should 
drop the assumption of homogeneity i.e., inhomogeneous solutions are useful. 

Szekeres [2] in 1975 gave a class of inhomogeneous solutions representing irrotational dust 
for the metric of the form (known as Szekeres metric) 



Subsequently, the solutions have been extended by Szafron [3] and Szafron and Wainwright 
[4] for perfect fluid and they studied asymptotic behaviour for different choice of the 
parameters involved. Later Barrow and Stein-Schabes [5] gave solutions for dust model 
with a cosmological constant and showed the validity of the Cosmic 'no-hair' Conjecture. 

In this work, we find inhomogeneous solutions for (n-|- 2)-dimensional Szekeres space-time 
with perfect fluid and a cosmological constant. The paper is organized as follows: The 
basic equations are presented in section II while the solutions have been written in section 
III. An asymptotic study of particular solutions are given in section IV. Finally the paper 
ends with a short discussion in section V. 



I. INTRODUCTION 



e^'-dr^ - e^^idx^ + dy^) 



II. BASIC EQUATIONS 



The metric for the (n + 2)-dimensional Szekeres space-time is in the form 



n 




(1) 



i=l 



where a and (3 are functions of all the [n + 2) space-time variables i.e., 



a = a(i, r, xi, ...,a;„), (3 = P{t,r,xi, ...,a;„). 
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The Einstein equations for the perfect fluid with a cosmological constant is of the form 

G^y = Agab + {p + p)uij.Uu - PQiiv (2) 

where p and p are energy density and isotropic pressure measured by an observer moving 
with with the fluid, A is the cosmological constant and is the fluid flow four vector. 
Since u = d/dt, the flow lines arc geodesies and the contracted Bianchi identities imply 
that pressure is a function of t only i.e., p = p{t). As there is no restriction on the energy 
density so p is in general a function of all the (n + 2) variables i.e., p = p{t, r, xi, a;„) and 
hence no equation of state is imposed. 

Now from the non-vanishing components of the field equations (2) for the above metric 
(1), we have 

1 " f 1 

nap + -n(n - l)/?^ - e'^/^ X] | < + 3 ('^ " " + - 2)a,,/3,, + a,,,, 

i=l 



+{n- } + e-^" <^ na'p' - -n{n+ l)/3'^ - n/3" }=A + p (3) 



ln{n + l)p^ + nP- ^n{n - l)e-^^P'^ - e-^^ ^ | i (n - 1) (n - 2)l3l^ + (n - | 

= K-p (4) 

a^+a + {n- l)d/3 + ^n{n - 1)/?^ + (n - + e'^" |(n - I) a! (3' - ^n{n - 1)/?'^- 
(n - 1)/?"} - ^ L^^ + l(n _ 2)(n - a)/^^^. + a,^,^ + (n - 2)/?,^,^. + (n - 3)a,^/3,^ j 

- e-^^ I (n - l)a,,/3,, + l(n - l)(n - 2)01 1 = A - p (5) 

Q!x3 (-Q!xi + /3xJ + /3x,- ("xi + (n - 2)/3xJ - a^ix^ - (n - 2)/3a:,^^. =0, (17^ j) (6) 

d/3' -jSP' -j3' =0 (7) 

- da^, + - d^. - (n - = (8) 



«x./3'-/3;. =0 (9) 

where dot, dash and subscript stands for partial differentiation with respect to t, r and 
the corresponding variables respectively (e.g., P^i = §§:) with i,j = 1,2, ...,n. 

Prom equations (7) and (9) after differentiating with respect to Xi and t respectively, we 
have the integrability condition 

/3xJ3" = 0, i = l,2,...,n (10) 

Thus, if /?' ^ we must have fixi = Oj « = 1,2, In the following section we shall 

consider the following possibilities 

(i)/3V0, (ii) /3' = 0, /j,, = 0, i = l,2,...,n 

to get solutions of the field equations. 



III. SOLUTIONS TO THE FIELD EQUATIONS 

In this section we shall solve the field equations (3)- (9) using the above restrictions 
separately. 

Case I: /?' ^ 

Here due to the restrictions (i^i =0, i = 1,2, ...,n we have from the field equations (7) 
and (9), the form of the metric coeSicient as 

e'^ = i?(f,r) e'^^'^'^i'-'^") (11) 

and 

e" = B! + Rv' (12) 

Now substituting these forms for the metric coefficient in equation (4) we have the differ- 
ential equations for R and v as 

RR+\{n- l)R? + ^{p{t) - A)R' = ^/(O (13) 

and 



^'"^ E { ('^ - + 2i^-.-i } = fir)-n (14) 



with /(r) as arbitrary function of r alone. 

Equation (13) can be integrated once to have the first integral 



where F{r) is another arbitrary function of r (appears due to integration). 

Also from (14) the solution for i/ will be 



n n 



e-- = A{r) E ^? + E ^» + (16) 



i=l i=l 

with the restriction 

n 

^i?2-4AC = /(r)-l (17) 

for the arbitrary functions ^(r), Bi{r), z = 1, 2, n and C(r). 
Further, to solve R completely let us consider p as a polynomial in t as 

p{t) = Pot-'' (18) 
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(po and a are positive constants) and we have the general solution for R as 



(19) 



according as ^ is an integer, non-integer and a = 2. Here Ci and C2 are arbitrary functions 
of r and we have chosen 



1 



a- 2' 



(n-l-l)po 1,, , n ^^ 



It is to be noted that to derive the above sohition wc have chosen A = = /(r). However, 
for non-zero A (but /(r) = 0) the solution is possible only for a = and 2 as 



' C\Cos{tJ^{po-m + C2Sin{tJ^{po-A)}, when a = 
R^ = \ (20) 

.V^{^^i^c[-^0^]+C2lc[-^V^l}' ^hena = 2 

with c = yi-'-^ . 

Further, if we consider the dust model (i.e.. p(t) = 0) then the above solution (19) 

simplifies to oc t. Hence for the usual 4D (i.e., n = 2) the scale factor R grows as 

^2/3 

as in the usual Priedmann model. 
Now, the physical and kinematical parameters have the following expressions 



_n F' + {n + l)Fu' _ po 
^~ 2 R"{R' + Ru') ~ ¥ 



(21) 



e = 



RR! + (n + l)RRv' + nRR' 
R{R' + Ru') 



(22) 



2E"-i(i?/' - 2i?7) + {n- 1){RF' - {n + 1)R'F) 



RR'^iR' + Rv') 



8(n+ l)(n- 1)2 
Case II: f3' = /j^. = 0, i = 1, 2, n 
In this case from the field equations we have the form of the metric functions 



(23) 



e'^ = R{t) 6-^(^1'^^'-'^") 



and 



(24) 



e" =R{t) r]{r,xi,X2,...,Xn) + IJ.{t,r) 



(25) 



Then as before from the field equation (4) we have similar differential equations in R and 
v as 
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and 



-1 -1 T7 — 1 

RR+-{n- 1)R' + -{p{t) - A)R^ = -^K 



e-^''Y.{{n-2)vl+2v,,,,}=K- 

1=1 



with K, an arbitrary constant. 

Here we take the solution for v in the form 



(26) 



(27) 



(28) 



where the arbitrary constants P,Qi (i = 1, 2, n) and S are restricted as before 



(29) 



Now to determine the function -q, we have from the field equation (6) 



52(e-^7y) 







dxidxi 

and then from the field equation (5) we have the solution 



(30) 



e "t] = u{r) a'i + Vi{r)xi + w{r) 



1=1 



with u{r),Vi{r) {i = 1,2, ...,n) and w{r) as arbitrary functions. 



(31) 



Further, to obtain the function ^ we use the following combination of the field equations 
namely. 



^G^:-G; = (n-l)(A-p(i)) 



i=l 



and the resulting differential equation in /x is 



Rfl + (n — l)Rfi + II 



R+- (p{t) - A)R 
n 



9{r) 



with 



(32) 



(33) 



g{r) = (n- 1) 



2{uS + wP)-^'> 



(34) 



For explicit solution if we choose p{t) as in the previous case (see eq.(18)) then the explicit 
form for R is same as in equation (19) except here Ci and C2 are arbitrary constants and 
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we have chosen K = A = 0. Similarly, we have the same solutions (20) for non-zero A. But 

wc note that the differential equation (33) is not solvable for any value of n. In fact only for 
n = 3 (i.e., for five dimension) we have the complete solution 



for a = 2 and ^2 = 5 ^1 + ^1 - with A = g{r) = 0. 

Further, the physical and kinematical parameters have the expressions as 



(35) 



2nfiR - K 
2/Ir 



n- 1 



T?(i?^/i - niiRR) - (n - l)ii'^R 
|JiR{^J. + riR) 



■P{t) 



(36) 



Rfi + ufiR + {n + l)RRri 
R{H + riR) 



(37) 



2(n + l) 



Rfi — R^ 
R{H + riR) 



(38) 



Finally, for simple dust case we have the solution for A = in terms of hypergeometric 
function as 



7^/, , X _ pi+i 4 n+1 3n-l (n + lf ^ 
Cx[t-to)-R^ 2ti[-,— — -,— — -, -r-7^J{r)R J 



^2'2n-2'2n-2' AnCi 
But for non-zero A we can have solution only for n = 3 as 



(39) 



i?2 = (716^^+^26- 



A 



However, we can have an integral equation from (15) (with p = 0) as 



(40) 



dR 



and we have the following particular solutions : 

(i) F(r) = 



(41) 



R- 



2A 



2A 



n{n+l) 



ill) f{r) = (i.e., if = in Case II) 



R^ 



{n + l)F(r) 
2A 



Sink 



{t - to) 



{n + l)A 
2n 



(iii) A = (n = 3) 
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(iv) /(r)=A = 

(v) /(r) = F(r) = 

^ ^ g(*-*o)\/^STT 

Also for the dust ji has the solution 

(42) 
for n = 3. 

IV. ASYMPTOTIC BEHAVIOUR 

We shall now discuss the asymptotic behaviour of the solutions presented in the previous 
section for both perfect fluid and dust model separately. The co-ordinates vary over the 
range : to <t <oo; — oo < r < oo; — oo < aij < oo, i = 1, 2, n. 

A. Perfect fluid model 

As p > 0, p ^ so we must have | < gi < 1. We shall first consider the case when a = 2. 
For large t {Case I i.e., /?' ^ 0) 

|(F' + (n + l)Fi/')i"'^^ 

For Case II, we have similar behaviour for large t together with /i ^ f^^ where qi is the 
value of qi for n = 3 and also we have ^ < gi < 1. Thus as i ^ oo, {p,p) fall off faster 
compare to {0,a), while the scale factor R (and fi) gradually increases with time. So the 
model approaches isotropy along fluid world line as t ^ oo. 



R = IJ.= 



+ C2e- 



K 
A 
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B. Dust model 

For the dust case with non-zero A (and n = 3) we have for large t, 

p = po, a constant 
= 9o, a constant 
a = ao, a constant 

We note that for Case I, po — ao — while for Case II we have po 7^ 0, ctq ^ 0. Thus 
universe will behave locally like de-sitter model in Case I though the global geometry 
will be different. However, in Case II the universe will not isotropize and for large t the 
measure of anisotropy and expansion scalar become finite constant. 



V. CONCLUDING REMARKS 

In this paper, we have found cosmological solutions for (n + 2)-dimensional Szekeres 
form of metric with perfect fluid (or dust) as the matter distribution. Wc can classify the 
solutions in two categories namely, (i) /3' ^ and (ii) /3' = 0. The first set of solutions are 
known as quasi-spherical solution while second class of solutions are termed as cylindrical 
type of solutions. 

However, if we assume the arbitrary functions A{r), Bi{r) and C(r) to have the constant 
values namely A{r) = C(r) = | and Bi{r) = (for all i = 1,2, n) i.e., if we have chosen 
the arbitrary function /(r) to be zero then using the transformation 

Xi = SinOnSinOn-i Sin92Cot^9\ 

X2 = Cos9nSin6n-i Sin92Cot^9i 

X3 = Cos9n-iSin9n-2 Sin62Cot^6i 



Xn-i = Cos63Sin92Cot^9i 
Xn = Cos92Cot^9\ 

the Szekeres metric (1) with the solution (11), (12), (16) reduces to the (n -|- 2)-dimensional 
spherically symmetric metric 

ds^ = - R'^dr^ - R^dnl 

It is to be noted that if the arbitrary functions A,Bi,C depend on r, then we can not get 

the above spherical form of the space-time. So these arbitrary functions play an important 
roll to identify the nature of the space-time. 

Finally, the study of asymptotic behaviour shows that some of the solutions will become 
isotropic at late time while there arc solutions for which shear scalar remains constant 
throughout the evolution. Thus Cosmic 'no-hair' Conjecture is not valid for all solutions. 
This violation of Cosmic 'no-hair' Conjecture is not unusual because the Szekeres metric 
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may not have always non-positive 3-space curvature scalar. Also in this higher dimensional 
Szckcrcs space-time we have solutions which expands as de Sitter in some directions but 
not in other directions. 
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